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Using  a  combination  of  electronic  -structure  theory,  variational  transition-state  theory,  and  solutions  to 
the  time-dependent  master  equation,  we  have  studied  the  kinetics  of  the  reaction  between  ethyl  and 
molecular  oxygen  theoretically  over  wide  ranges  of  temperature  and  pressure.  The  agreement  between 
theory  and  experiment  is  quite  good.  By  comparing  the  theoretical  and  experimental  results  describing 
the  kinetic  behavior,  we  have  been  able  to  deduce  a  value  for  the  C2H5-02  bond  energy  of  —34  kcal/mol 
and  a  value  for  the  exit-channel  transition-state  energy  of  —4.3  kcal/mol  (measured  from  reactants).  These 
numbers  compare  favorably  with  our  electronic-structure  theory  predictions  of  33.9  kcal/mol  and  —3.0 
kcal/mole,  respectively.  The  master-equation  solutions  show  three  distinct  temperature  regimes  for  the 
reaction,  discussed  extensively  in  this  paper.  Above  T  ««  700  K,  the  reaction  can  be  written  as  an  elementary 
step,  C2H5  +  02  <-»  C2H4  +  H02,  with  the  rate  coefficient  k(T)  =  3.19  X  10~17  T102  exp(2035/RT) 
cm3  (molecules  s)  independent  of  pressure,  even  though  the  intermediate  collision  complex  may  suffer  a 
large  number  of  collisions. 


Introduction 

The  reactions  of  alkyl  radicals  with  molecular  ox¬ 
ygen  constitute  an  important  class  of  reactions  in 
combustion.  It  is  the  properties  of  these  reactions 
that  cause  the  change  in  oxidation  mechanism  of  al¬ 
kanes  between  500  K  and  800  K  and  the  “negative 
temperature  coefficient”  of  their  overall  oxidation 
rate  in  the  same  temperature  regime  [1—3].  Ethyl 
(C2H5)  is  the  prototype  alkyl  radical,  and  conse¬ 
quently  its  reaction  with  02  has  been  studied  exten¬ 
sively  in  the  laboratory  [4-12],  There  have  also  been 
several  attempts  to  treat  the  reaction  theoretically 
[8,13,14],  the  most  extensive  of  which  is  the  study 
of  Wagner  et  al.  [8],  The  latter  authors  provided  a 
comprehensive  review  of  the  literature  through 
1990,  so  we  may  dispense  with  that  here. 

The  present  investigation  is  theoretical.  It  is  a 
rather  substantial  extension  of  the  work  of  Wagner 
et  al.,  who  treated  the  C2H5  +  02  reaction  as  the 
sum  of  four  separate  elementary  reactions,  each  with 
its  own  pressure-  and  temperature-dependent  rate 
coefficients.  These  rate  coefficients  were  calculated 
from  chemically  activated  Rice-Ramsperger-Kas- 
sel-Marcus  (RRKM)  theory  using  canonical  varia¬ 
tional  transition-state  theory  for  the  loose  transition 
state  and  a  pseudo  strong-collider  model  to  calculate 


the  effects  of  collisions.  In  the  present  investigation, 
we  use  ab  initio  electronic  structure  theory  (density 
functional  theory  [DFT]  and  a  G2-like  method)  to 
characterize  various  features  of  the  potential  energy 
surface,  and  we  use  microcanonical  (/iVT)  and  mi- 
crocanonical/J-conservative  (/iVT-J)  variational  tran¬ 
sition-state  (RRKM)  theory  to  calculate  the  high- 
pressure  limit  of  the  addition  reaction  and  the 
zero-pressure  limit  of  the  bimolecular  channel.  Most 
importantly,  however,  we  use  solutions  to  the  time- 
dependent  master  equation  to  characterize  the  re¬ 
action  over  wide  ranges  of  temperature  and  pres¬ 
sure,  including  the  “transition  region”  between  500 
K  and  800  K. 


Theory 

Figure  1  shows  a  potential  energy  diagram  for  the 
C2H5  +  02  reaction.  The  energies  of  the  stationary 
points  on  the  potential  energy  surface  (PES)  come 
from  the  work  described  below.  Our  objective  is  to 
calculate  the  total  rate  coefficient  for  the  reaction 
k(T,p),  and  the  branching  fraction  a{T,p),  that  is, 
the  fraction  that  forms  the  bimolecular  products 

c2h4  +  ho2. 
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Fig.  1.  Reaction  coordinate  diagram  for  the  C2H5  +  02 
reaction,  based  primarily  on  the  G2-like  energies  of  Table 
1.  The  energies  shown  for  the  stationary  points  of  the  po¬ 
tential  energy  surface  are  potential  energy  plus  zero-point 
vibrational  energy. 


TARLE  1 

Energetics  relative  to  reactants  for  the  stationary 
points  in  the  C2H5  +  02  reaction  (energies  in  keal/ mol) 


Species 

B3LYP 

G2c-like 

<S2>d 

exp® 

C2H502(g) 

-29.0“ 

-33.9 

0.763 

-34.0 

C2H502(t) 

-29.0“ 

-33.8 

0.763 

-34.0 

TS-2 

-1.9“ 

-3.0 

0.855 

-4.3 

TS-3 

8.0“ 

3.1 

0.813 

CH2CH2OOH 

-9.9“ 

- 17.0 

0.763 

TS-5 

1.5“ 

1.9 

1.023 

CH2CH2  •  •  •  HOO 

- 14.0“ 

-16.6 

0.762 

c2h4  +  ho2 

-11.2“ 

-12.7 

0.763 

TS-6 

1.5b 

-0.6 

1.297 

C2H40  +  OH 

—  25.2b 

-34.0 

0.757 

TS-4 

14.8b 

8.2 

0.950 

CH3CHO  +  OH 

-  56.5b 

-61.3 

0.757 

Abstraction 

9.5b 

18.2 

1.687 

“From  Ref.  [16]  using  the  TZ2Pf  basis  function. 

bPresent  work  employing  the  6-311  +  +  G(d,  p)  basis  set. 

“Approximate  E[QCISD(T)/6-311  +  +G(3df  2 pd)]  en¬ 
ergy  as  described  in  equation  1. 

Expectation  value  of  the  total  electron  spin  squared  for 
the  MP2  wave  function  used  in  the  G2-like  calculations. 

“Deduced  from  comparing  predicted  kinetic  behavior 
with  experiment. 

Quantum  Chemistry 

The  C2H5  +  02  reaction  has  been  the  subject  of 
a  number  of  ab  initio  quantum  chemical  studies  (cf. 
[15-17]  and  references  cited  therein)  with  methods 
ranging  from  density  functional  theory  to  more  tra¬ 
ditional  MPn,  coupled  cluster,  and  configuration  in¬ 
teraction  methods.  The  vibrational  frequencies  ob¬ 
tained  in  the  B3LYP  [18J/DZP  analysis  of  Ref.  [16] 


are  employed  in  the  present  investigation.  In  the  in¬ 
terest  of  further  delineating  the  allowable  energet¬ 
ics,  we  have  evaluated  stationary-point  electronic  en¬ 
ergies  with  an  approach  that  is  closely  analogous  to 
the  G2  method  [19]: 

E[QCISD(T)/6-311  +  +G(3 df  2 pd)] 

=  £[QCISD(D/6-311+  +  G(d,p)] 

+  E[MP2/6-311+ +G(3 df,  2 pd)] 

-  £[MP2/6-311+  +G(d,p)] 

+  HLC  +  E(0)  (1) 

where  E®  is  the  zero-point  vibrational  energy. 

These  evaluations  and  the  rate  coefficient  calcu¬ 
lations  discussed  below  employ  molecular  structures 
obtained  from  B3LYP/6-311  +  +G(d,p)  optimiza¬ 
tions.  The  higher-level  corrections  (HLCs)  used  in 
different  variants  of  G2  theory  [19-21]  yield  a  de¬ 
crease  in  the  energy  (relative  to  reactants)  of  about 
2.6  to  3.1  kcal/mol  for  each  of  the  stationary  points 
of  interest  here.  In  contrast,  the  more  recent  G3 
theory  [22]  suggests  the  HLCs  should  be  only  —  0.3 
kcal/mol,  and  that  this  is  likely  to  be  a  better  cor¬ 
rection  even  for  calculations  at  the  G2  level.  Thus, 
this  smaller  correction  is  employed  here. 

The  results  of  these  G2-like  calculations  are  pro¬ 
vided  in  Table  1  and  are  seen  to  be  in  reasonable 
agreement  with  the  corresponding  B.3LYP  calcula¬ 
tions  of  Ref.  [16],  but  they  do  indicate  a  general 
lowering  of  the  energies.  Also  shown  in  Table  1  are 
the  stationary-point  energies  deduced  by  comparing 
our  rate-coefficient  calculations  with  experimental 
results  (discussed  below).  These  energies  are  in  good 
agreement  with  the  ab  initio  results. 

The  Loose  Transition  State 

There  are  no  electronic  structure  calculations 
from  which  to  draw  information  about  the  PES 
along  the  C2H5  +  02  entrance  channel,  and  it  is 
beyond  the  scope  of  the  present  work  to  provide 
such  information.  Consequently,  we  approximated 
the  potential  in  this  region  of  the  surface  by  the 
method  described  by  Miller  and  Klippenstein  [23], 
Briefly,  the  potential  consists  of  three  parts. 

1.  The  potential  along  the  reaction  coordinate. 
This  part  of  the  potential  is  approximated  by  the  Var- 
shni  function  [23,24],  The  two  parameters  in  the 
Varshni  potential  come  from  knowledge  of  the  bond 
energy  and  from  a  DFT  calculation  of  the  force  con¬ 
stant  matrix  at  the  potential  minimum. 

2.  The  potential  for  the  “conserved”  degrees  of 
freedom  orthogonal  to  the  reaction  coordinate.  This 
part  of  the  potential  corresponds  to  the  degrees  of 
freedom  that  can  be  identified  as  normal-mode  vi¬ 
brations  in  the  separated  fragments  and  is  assumed 
to  be  the  same  as  in  the  fragments. 
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3.  The  potential  for  the  “transitional”  degrees  of 
freedom  orthogonal  to  the  reaction  coordinate.  This 
part  of  the  potential  is  described  in  terms  of  a  set  of 
internal  angles.  The  potential  is  written  essentially 
as  a  sum  of  products  (in  pairs)  of  sinusoidal  func¬ 
tions,  with  the  appropriate  phases  and  periods  de¬ 
termined  by  symmetry  (see  equation  4  of  ref.  23). 
The  coefficients  in  the  expression  are  functions  of 
the  bond  distance  R  and  are  determined  by  calcu¬ 
lating  the  appropriate  force  constant  matrix  using 
DFT  at  the  potential  minimum,  Fy(R0)  and  assum¬ 
ing  that  these  matrix  elements  decay  exponentially 
with  bond  distance: 

FfR)  =  FfR0)  exp [~r,(R  -  R0)]  (2) 

The  only  adjustable  constant  is  //,  and  we  adjust  its 
value,  rj  =  1.45  A-1,  to  give  an  accurate  represen¬ 
tation  of  the  high-pressure  (p),  low-temperature  ( T ) 
rate  coefficients. 

The  sum-of-states  functions  Nf(EJ)  and  JVf  (£) 
are  calculated  with  the  VARIFLEX  [25]  computer 
code  using  methods  developed  by  Klippenstein  [26]. 
In  the  former  case,  there  is  one  transition-state  di¬ 
viding  surface,  TS-1,  for  every  EJ  combination,  and 
in  the  latter  there  is  one  for  every  E,  where  E  is  the 
total  energy  and  J  is  the  total  angular  momentum 
quantum  number. 


Hindered  Rotors,  Symmetry  Numbers,  and  State 
Counting 

The  CH3  ■  ■  ■  CH2  torsional  mode  in  the  ethyl  rad¬ 
ical  is  treated  in  the  present  investigation  as  a  hin¬ 
dered  rotor  using  a  Pitzer-Gwinn  approximation  [27] 
for  the  canonical  partition  function  of  the  reactants. 
The  corresponding  contribution  to  the  transition- 
state  partition  function  has  been  evaluated  with  a 
microcanonical  analog.  In  particular,  the  hindered- 
rotor  density  of  states  was  evaluated  as  the  classical 
hindered  rotor  density  times  the  quantum  harmonic 
density  divided  by  the  classical  harmonic  density. 
The  hindered-rotor  rotational  constant  (15.2  cm-  ) 
was  evaluated  according  to  the  prescription  provided 
by  Pitzer  in  equation  1  of  Ref.  [28].  A  sinusoidally 
hindered  potential  of  period  2n/3  was  employed 
with  the  barrier  height  set  to  yield  the  ah  initio  har¬ 
monic  vibrational  frequency  at  the  minimum  of  the 
hindering  potential.  Rotational  symmetry  numbers 
of  3  and  2  were  employed  for  the  ethyl  and  oxygen 
fragments,  respectively.  The  transition-state  evalua¬ 
tions  for  the  entrance  channel  included  the  contri¬ 
bution  from  the  binding  of  only  one  of  the  two  ox¬ 
ygens.  Thus,  the  net  symmetry  number  for  the 
entrance  channel  transition  state  was  reduced  from 
6  to  3.  The  reaction  was  assumed  to  occur  only  on 
the  ground  doublet  electronic  PES,  and  so  an  elec¬ 
tronic  correction  factor  of  1/3  was  applied  to  the 
total  rate  coefficient. 


For  the  complex  (CH3CH202),  we  have  treated 
both  the  OO  ■  •  •  CH2CH3  and  the  CH3  •  •  •  CH200 
torsional  modes  as  hindered  rotors  using  a  micro- 
canonical  analog  of  the  Pitzer-Gwinn  approximation. 
The  internal  rotational  constants  for  these  two  hin¬ 
dered  rotors  are  2.22  cm-1  and  5.65  cm-1,  respec¬ 
tively.  Hindering  potentials  of  period  2n/3  were  em¬ 
ployed  for  both  these  modes,  since  the  trans  and  two 
gauche  configurations  are  accidentally  degenerate 
(to  within  0.1  kcal/mol).  The  combination  of  this  ac¬ 
cidental  degeneracy  and  the  enantiomeric  relation 
of  the  two  gauche  configurations  yields  an  effective 
rotational  symmetry  number  of  3.  The  density  of 
states  of  the  complex  pc(E,J)  was  calculated  by  con¬ 
volving  the  exact-count  harmonic  contributions  with 
both  the  internal  and  external  rotations.  When 
needed,  pc(E)  was  calculated  as  (2 J  +  1  )pc(E,J). 

For  the  exit  channel  transition  states  (TS-2,  TS-3, 
and  TS-4),  all  vibrational  modes  were  treated  as  har¬ 
monic  oscillators,  and  the  rotational  symmetry  num¬ 
ber  was  unity.  The  sums  of  states  Nf(E,J),  Nf(E,J), 
and  Nf(E,J)  were  computed  from  conventional 
transition-state  theory  by  the  usual  methods  [29], 
The  corresponding  N±(E)  values  were  calculated 
from  N±(E)  =  (2/  +  1  )N±(E,J). 


Rate  Coefficient  Calculations 

First,  consider  the  two  limits,  the  high-pressure 
limit  of  the  addition  reaction  and  the  zero-pressure 
limit,  which  of  course  includes  only  the  bimolecular 
product  channel.  Using/; VT-J,  we  can  write  the  cor¬ 
responding  rate  coefficient  expressions  as 


X  exp  ( —  fE)dE  (3) 

and 


X 


.3 


Nf(E,J)Nj(E,J) 
Nf(E,J)  +  Nf(E,J) 


exp(— ySE)dE 


(4) 


where  f  =  (kBT)_1,  kB  is  Boltzmanns  constant,  h  is 
Planck’s  constant,  and  ()r  is  the  reactant  partition 
function,  including  relative  translational  contribu¬ 
tions.  In  these  expressions,  the  electronic  degener¬ 
acies  and  symmetry  numbers  have  been  subsumed 
into  <?r  and  the  N±(E,J)  s.  The  /<VT  approximations 
to  k0(T)  and  k„(T)  are  similar  to  the  above  except 
that  there  is  no  sum  over  J,  the  N±’s  are  only  func¬ 
tions  of  E,  and  there  is  only  one  TS-1  for  eveiy  E. 
Equation  4  assumes  that  TS-2  gives  the  only  contri¬ 
bution  to  the  bimolecular  channel  and  that,  once  this 
TS  is  traversed,  a  complex  inevitably  goes  on  to  form 
C2H4  +  H02.  The  latter  is  a  safe  assumption,  and 
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Fig.  2.  Rate  coefficients  as  a  function  of  temperature  for 
various  pressures. 


to  check  the  former  we  have  performed  a  few  cal¬ 
culations  with  N%(EJ)  in  equation  4  replaced  by 
N£(EJ)  +  Nf(EJ)  +  N£(EJ).  This  gives  us  an 
estimate  of  the  upper  limit  of  the  error  incurred  by 
neglecting  TS-3  and  TS-4  as  possible  exit  paths  from 
the  initial  complex. 

To  calculate  rate  coefficients  as  a  function  of  T  and 
p,  we  solve  the  time-dependent,  one-dimensional 
master  equation.  To  formulate  the  problem,  envision 
the  situation  typical  in  kinetics  experiments, 
«c2h5  «o2  <  »He  where  n  denotes  number  den¬ 

sity.  For  the  present  investigation  we  assumed  that 
the  bath  gas  is  always  helium  and  that  the  reactants 
are  constantly  maintained  in  a  thermal  distribution. 
Under  these  conditions  we  can  write  a  linear  master 
equation  for  n{E),  where  n(E)dE  is  the  number  den¬ 
sity  of  C2H502  complexes  between  energies  E  and 
E  +  dE,  and  a  linear  rate  equation  for  «c2h5-  Space 
limitations  preclude  giving  a  derivation  here,  but  it 
is  possible  to  express  these  equations  in  a  simple 
form.  After  approximating  the  collision  integral  in 
the  master  equation  as  a  sum,  the  problem  can  be 
posed  in  Dirac  notation  as 

\y(t))  =  G\ij(t))  (5) 

In  equation  5,  I ij(t))  is  the  vector  of  unknowns,  I (/(f)) 

— >  Vji,  •  •  •  ,  iji,  ■  ■  •  ,  iJm,  ( Keqii0J d£)1/2xC2H5]  , 

where  iji  =  ij(Ej ),  Et  =  E0  +  (i  —  I)  AE,  i  =  1,  . 

M,  M  is  the  number  of  grid  points  in  the  energy 
space,  AE  is  the  spacing  between  grid  points,  E0  is 
the  energy  corresponding  to  the  ground  vibrational 
state  of  C,H502,  xC2h5  =  «c2h5(1)/«c2h5(0), 
ij(E)  =  x(E)/VF(E),  x(E)  =  n(£)/»C2H  (0),  F(E)  is 
the  equilibrium  energy  distribution  of  the  C2H502 
adduct,  and  is  the  equilibrium  constant  for  the 
addition  reaction,  C2H5  +  02  ++  C2H502.  The  op¬ 
erator  G  in  the  discrete  energy  basis  is  a  real  sym¬ 
metric  (Hermitian)  matrix,  which  allows  a  simple  so¬ 
lution. 


One  may  write  the  solution  to  equation  5  as 

M+l 

l;/(f)>  =  2  eVlgi>  <gili/(0)>  (6) 

i  =  1 

where  and  Igj)  are  the  M+l  (negative)  eigen¬ 
values  and  eigenvectors  of  G,  and  li/(0)}  corresponds 
to  xC2Hs(0)  =  1.  Methods  for  diagonalizing  real, 
symmetric  matrices  are  highly  developed,  and  good 
software  is  relatively  easy  to  find.  We  used  the 
DSYEV  routine  from  LAPACK  [31]  to  diagonalize 
G. 

From  the  solution  vector  \y(t)),  one  obtains 
xc2h5  (t)  —  [K-e^noJAE)  1/2  ijm+i  (t) 

Xj(t)  =  Jf{Ej)  y/f),  j  =  1, ....  M  (7) 

M 

*c2h5o2  (t)  =  2  Xi(t)AE 
1  =  1 

and 

*P  (t)  ~  1  —  xc2h5  if)  ~  xc2h5o2  W 

where 


_  kC2Hj  it)  _  »HQ2  it) 

F  nC2H5  (0)  «C2H5  (0) 

In  interpreting  our  results,  we  calculate  k(T,p,t), 
which  is  defined  as 


k(T,  p ,  t) 


1  <frc2H5  it) 

no2xc2H5  (1)  dt 


(8) 


If  all  the  C2H5  ultimately  reacts  (it  does  in  the  cases 
considered  here),  and  there  are  good  exponential 
decays,  k(T,p,t)  —  k(T,p)  (the  total  rate  coefficient) 
is  constant  in  time.  At  sufficiently  high  pressure,  in 
the  temperature  range  roughly  between  575  K  and 
700  K,  multiple  exponential  decays  occur.  In  these 
cases  we  characterize  the  rate  using  equation  8  at 
the  time  when  Xc2u5  =  0.5.  When  there  is  any  am¬ 
biguity  in  the  branching  fraction,  we  define  a(T,p), 
the  fraction  of  the  reaction  that  produces  C2H4  + 
H02,  as  Xp(f)  at  a  time  when  xC2u5  *=  0.01. 

In  all  our  master  equation  calculations,  we  used  a 
standard  exponential-down  model  for  the  energy 
transfer  function  with  ( AE j)  =  200  cm-1.  All  the 
calculations  were  done  with  VARIFLEX. 


Results  and  Discussion 

First,  we  consider  the  behavior  of  the  reaction 
from  a  broad  perspective.  In  Fig.  2,  k(T)  is  plotted 
for  various  pressures,  along  with  a  number  of  exper¬ 
imental  results  [6,8,10,11]  at  pressures  ranging  from 
slightly  less  than  1  torr  to  almost  an  atmosphere.  The 
upper  and  lower  solid  curves  represent  k^T)  and 
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Time  (milliseconds) 

Fig.  3.  Plot  of  xC2H5(f),  xC2h5o2(°°)  “  xc2h5o2W,  ancl 
Xp(°°)  —  xp (t)  on  a  semi-log  plot  for  T  =  300  K,  p  =  1 
torr.  Straight  lines  indicate  “normal”  behavior. 


Fig.  4.  Comparison  of  high-pressure  rate  coefficients 
with  experiment. 


k0(T)  respectively,  calculated  from  pVT-J.  As  ex¬ 
pected,  except  for  a  pathological  calculation  at  700 
K  and  10  atm  (where  it  was  necessary  to  use  the 
somewhat  arbitrary  definition  of  k(T,p)  discussed 
above),  till  the  rate  coefficients  lie  between  these  two 
boundaries  (within  reasonable  uncertainty  for  the 
experiments).  From  the  master  equation  calcula¬ 
tions  it  is  possible  to  identify  three  different  regimes 
of  the  reaction. 

1.  Low-temperature  regime.  In  this  regime  the  re¬ 
action  exhibits  normal  behavior.  The  rate  coefficient 
and  branching  fraction  are  functions  of  T  and  p,  and 
the  time  evolution  is  perfectly  exponential  (i.e.,  the 
rate  coefficient  is  well  defined).  This  regime  persists 
up  to  T  ~  600  K. 

2.  Transition  regime.  Roughly  between  600  K  and 
700  K  the  xC2H5(f)  is  decidedly  non-exponential,  at 
least  at  sufficiently  high  pressure.  Consequently,  it  is 
impossible  to  extract  a  rate  coefficient  rigorously 
from  the  calculation,  at  least  in  the  normal  sense. 
The  temperatures  at  which  this  occurs  correspond 
to  the  equilibration  of  the  hypothetical  C2H5  +  Oz 

C2H502  reaction.  For  the  conditions  of  our  cal¬ 
culations,  the  equilibrium  fraction  of  C2H5  tied  up 


as  C2H502  at  500  K  is  0.997,  but  it  drops  off  precip¬ 
itously  to  0.108  at  700  K  and  to  0.0051  at  800  K. 
The  complex  behavior  of  the  reaction  in  this  regime 
is  due  to  the  simultaneous  dissociation  of  C2H502  to 
both  C2H5  +  02  and  C2H4  +  H02  after  varying 
degrees  of  stabilization.  Note  the  dramatic  drop-off 
of  the  rate  coefficient  in  Fig.  2  at  high  pressures 
between  650  K  and  700  K  as  the  stabilization  reac¬ 
tion  equilibrates.  Stabilization  becomes  problematic 
under  these  conditions,  and  the  only  outlet  is  the 
bimolecular  channel.  The  temperature  range  for  this 
regime  is  sensitive  to  the  C2H5-02  bond  energy.  We 
chose  the  value  of  34  kcal/mol  used  in  the  calcula¬ 
tions  to  agree  with  experimental  observations.  It  is 
reassuring  that  our  G2-like  calculation  of  33.9  kcal/ 
mol  is  in  excellent  agreement  with  this  number.  Pre¬ 
vious  estimates  [8,30]  are  also  compatible  with  these 
results,  but  the  DFT  bond  energy  of  29  kcal/mol  is 
clearly  too  small. 

3.  High-temperature  regime.  The  onset  of  this  re¬ 
gime  is  at  T  =  700  K.  By  750  K,  the  only  products 
of  the  reaction  are  C2H4  +  H02,  and  the  rate-co- 
efficient  curves  k(T,p)  in  Fig.  2  for  all  pressures  have 
coalesced  to  the  zero-pressure  limit.  This  latter  re¬ 
sult  is  somewhat  surprising,  particularly  when  one 
realizes  that  the  C2H502  complex  still  suffers  nu¬ 
merous  collisions  at  high  pressure.  In  this  regime, 
the  reaction  can  be  characterized  by  a  single  set  of 
products,  C2H4  +  H02,  and  a  single  rate  coefficient, 
k(T),  independent  of  pressure.  As  the  last  sentence 
implies,  one  again  obtains  good  exponential  decays 
in  *c2hsM- 

Figure  3  demonstrates  what  we  mean  by  normal 
behavior.  It  is  a  plot  at  T  =  300  K,  p  =  1  torr  of 
xc2h5(/)>  xc2h5o2(C0)  —  xc2h5o2W>  and  xp(°°)  —  xp(f) 
on  a  log  scale.  If  the  reaction  is  governed  by  a  single 
time  constant  (i.e.,  a  single  eigenvector  of  G),  and 
one  can  thus  identify  a  “good”  rate  coefficient,  all 
these  curves  should  be  straight  lines.  They  are,  at 
least  on  the  time  scale  for  C2H5  removal. 

Figure  4  is  a  comparison  of  our  high-p,  low-T  rate 
coefficients  with  the  experimental  results  of  Kaiser 
and  coworkers  [9,11].  The  agreement  between  the¬ 
ory  and  experiment  is  excellent.  At  580  torr,  both  the 
experimental  and  theoretical  rate  coefficients  re¬ 
main  relatively  constant  from  250  K  up  to  T  ^  400 
K,  and  then  they  drop  off  more  rapidly  with  T.  This 
constancy  of  k(T ,  580  torr)  with  T,  even  with  a  weak- 
collider  bath  gas  such  as  helium,  is  a  strong  indica¬ 
tion  that  k„(T)  rises  with  temperature,  an  unusual 
result  (although  not  unheard  of)  for  a  barrierless  re¬ 
action.  Our  pVT-J  prediction  of  k„(T)  falls  within  the 
error  bars  of  the  room  temperature  determination 
of  this  rate  coefficient  reported  in  Ref.  [9],  Wagner 
et  al.  [8]  give  k„(T)  “7  X  10-12  cm3/(molecule-s) 
almost  independent  of  temperature,  a  value  some¬ 
what  smaller  than  the  experimental  results  and  the 
present  theoretical  predictions.  It  is  also  worth  not¬ 
ing  that  the  present  pVT-J  and  ju\ T  predictions  of 
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Fig.  5.  Comparison  of  low-pressure  rate  coefficients 
with  experiment. 


Fig.  6.  Comparison  of  the  branching  fraction  a(T)  with 
the  experiments  of  Clifford  et  al.  [12]  for  a  fixed  helium 
number  density  of  1.1  X  10ls/cm3. 


Fig.  7.  Time  history  of  xColI(t ),  xCilo„(t),  and  1  -  x  (t) 
for  T  =  700  K  ,p  =  100  torr. 


k^iT)  do  not  differ  significantly — the  /iVT  rate  co¬ 
efficient  is  only  6%  larger  at  250  K  and  impercep- 
tively  larger  for  T  >  400  K.  The  present  prediction 
for  fc„(T)  can  be  expressed  reasonably  accurately  in 
the  modified  Arrhenius  form,  kjT)  =  4.84  X  10” 13 
T0'32  cm3/(molecule-s)  for  250  K  <  T  £  2000  K. 

Figure  5  is  a  comparison  of  our  theoretical  pre¬ 
dictions  with  the  low-pressure  rate  coefficients  mea¬ 
sured  by  Gutman  and  collaborators  [6-8],  All  their 


experiments  were  conducted  at  pressures  ranging 
from  slightly  less  than  1  torr  to  approximately  15  torr. 
Our  theoretical  curves  corresponding  to  p  ~  I  torr 
and  p  —  20  torr  should  more  or  less  bracket  these 
results.  Within  reasonable  experimental  uncertainty, 
they  do.  The  rate  coefficients  in  the  high-tempera¬ 
ture  regime  (T  >  700  K)  are  governed  almost  exclu¬ 
sively  by  TS-2,  and  we  adjusted  the  value  of  its  en¬ 
ergy,  E2,  in  our  calculations  to  give  good  agreement 
with  the  experimental  rate  coefficients  in  this  re¬ 
gime.  This  value  of  E2  —  —  4.3  kcal/mol  (measured 
from  reactants)  is  in  reasonably  good  agreement  with 
our  G2-like  ab  initio  result  of  E2  =  —  3.0  kcal,  but 
it  differs  somewhat  from  the  DFT  calculation  of 
—  1.9  kcal/mol.  Using  E2  —  —  1.9  kcal/ mole  in  the 
calculations  gives  rate  coefficients  in  the  high-tem- 
perature  regime  that  are  substantially  smaller  than 
those  shown  in  the  figure  (by  approximately  a  factor 
of  4  at  T  =  1000  K). 

Also  shown  in  Fig.  5  are  the  fi\ T  and/iVT-J  results 
for  k0(T).  The  pVT  rate  coefficient  is  29%  larger  than 
that  calculated  from  p  VT-J  at  T  —  250  K,  a  differ¬ 
ence  that  diminishes  with  increased  temperature  un¬ 
til  it  is  virtually  imperceptible  for  T  >  500  K.  The 
two  rate  coefficients  come  together  at  high  tem¬ 
perature  because  k0(T)  begins  to  be  controlled  ex¬ 
clusively  by  the  single  tight  transition  state  TS-2,  that 
is  the  flux  through  TS-1  becomes  effectively  infinite 
(see  equation  4),  and  we  are  left  with  a  conventional 
transition-state  theory  result,  with  TS-2  the  transi¬ 
tion  state.  The  contribution  to  k0(T)  from  TS-3  and 
TS-4  is  no  more  than  9%  at  2000  K  (calculated  using 
the  G2-like  barrier  heights  of  Table  1);  it  is  much 
less  at  lower  temperatures.  For  700  K  <  T  <  2000 
K,  we  can  express  k0(T)  as  k0(T)  =  3.19  X  10”17 
T1 02  exp  (2035 /RT)  cm3/(molecule-s),  which  should 
be  sufficiently  accurate  for  flame  modeling. 

The  transition-state  energy  E2  also  plays  a  major 
role  in  determining  the  branching  fraction  cn(T,p)  in 
the  low-temperature  regime.  Fig.  6  compares  our 
predictions  for  a  with  the  experimental  results  of 
Clifford  et  al.  [12],  who  measured  the  branching 
fraction  from  room  temperature  to  700  K  at  a  con¬ 
stant  helium  number  density  of  1.1  X  1018/cm3.  The 
agreement  between  theory  and  experiment  is  quite 
satisfactoiy.  Both  show  a  slow  rise  in  a  from  room 
temperature  to  about  575  K,  followed  by  the  very 
rapid  rise  between  575  K  and  700  K  that  occurs  as 
the  C2H5  +  Oz  C2H502  reaction  equilibrates. 

Figure  7  shows  the  time  history  of  xCoH5(f), 
xc2h5o2(/);  and  1  —  xp (/)  on  a  semi-log  plot  for  T  = 
700  K,  p  =  100  torr.  This  plot  is  characteristic  of 
the  behavior  we  observe  throughout  the  transition 
regime  (at  sufficiently  high  pressure).  The  reactant 
concentration  shows  a  distinctly  non-exponential  de¬ 
cay.  The  function  .vC2h5(/)  first  decays  rapidly  as 
xc2h5o2W  rises.  After  xC2H  o2(f)  peaks,  both  xC2h5(/) 
and  Xc2H  o  Jf)  decay  with  the  same  exponential  time 
constant.  The  most  interesting  result  from  Fig.  7  is 
that  the  bimolecular  products,  xp(f),  are  governed  by 
this  same  long-term  time  constant  even  at  very  short 
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times.  However,  at  lower  pressures  xp(t)  can  also  dis¬ 
play  two-time-constant  behavior,  more  pronounced 
in  the  600  K  to  650  K  temperature  range  than  at  700 
K. 

It  is  instructive  to  interpret  our  results  in  terms  of 
the  eigenvalues  and  eigenvectors  of  G.  Two  eigen¬ 
values  play  a  major  role.  At  low  temperature,  it  is 
the  second  largest  eigenvalue  of  G  (i.e.,  the  second 
least  negative),  A2,  that  determines  the  time  evolu¬ 
tion  of  both  the  reactants  and  products.  The  largest 
eigenvalue,  A1;  is  associated  with  the  thermal  disso¬ 
ciation  of  C2H502  to  C2H4  +  HOz  and  is  much 
larger  than  A2.  This  largest  eigenvalue  decreases  rap¬ 
idly  with  increasing  temperature  (at  least  above  500 
K),  and  in  the  transition  regime  A4  and  A2  are  of 
comparable  size.  We  thus  obseive  the  behavior  de¬ 
scribed  above.  In  the  high-temperature  regime, 
there  is  again  a  separation  of  the  eigenvalues,  but 
now  it  is  A4  that  governs  the  time  evolution  of  C2Hs 
and  02  to  C2H4  +  H02. 

A  question  about  the  importance  of  the  direct  ab¬ 
straction  channel  always  arises  in  discussions  of  alkyl 
radical  +  02  reactions.  For  the  present  case,  we 
have  used  the  G2-like  transition  state  of  Table  1  to 
calculate  a  direct  abstraction  rate  coefficient 

kJJ)  =  1.93  X  10-17  T1'57  exp(  —  20578/RT) 

for  1000  K  <  T  £  2000  K.  This  rate  coefficient  is 
about  a  factor  of  10  smaller  at  2000  K  than  the  k(T) 
discussed  above.  The  difference  is  even  greater  at 
lower  temperatures,  and  thus  it  is  probably  reason¬ 
able  to  neglect  fca(T)  in  most  applications. 

Concluding  Remarks 

We  have  used  a  variety  of  theoretical  techniques 
to  study  the  reaction  between  ethyl  and  02  over  very 
wide  ranges  of  temperature  and  pressure.  These  in¬ 
clude  electronic-structure  theory,  variational  transi¬ 
tion-state  theory,  and  solutions  to  the  time-depen¬ 
dent  master  equation.  Our  predictions  of  rate 
coefficients  and  product  distributions  are  in  very 
good  agreement  with  experiment.  Moreover,  the  val¬ 
ues  of  the  C2H5-02  bond  energy  and  the  exit-chan¬ 
nel  transition-state  energy  of  34  kcal/mol  and  —4.3 
kcal/mol  (measured  from  reactants),  respectively, 
deduced  from  the  master  equation  results,  compare 
favorably  with  the  values  of  33.9  kcal/mol  and  —  3.0 
kcal/mol  calculated  from  our  G2-like  method. 

The  most  interesting  aspect  of  our  analysis  is  the 
prediction  of  three  different  regimes  of  the  reaction, 
discussed  extensively  in  the  text.  For  temperatures 
larger  than  T  700  K,  the  reaction  can  be  written 
as  a  simple  elementary  step,  C2H5  +  02  C2H4 
+  H02,  with  a  rate  coefficient  k(T)  independent  of 
pressure,  even  though  the  reaction  goes  through  an 
intermediate  complex  that  may  suffer  numerous  col¬ 
lisions.  Moreover,  the  rate  coefficient  is  the  same  as 


that  at  the  zero-pressure  (collisionless)  limit.  In 
modified  Arrhenius  form,  we  can  write  k(T)  in  this 
regime  as  k(T)  =  3.19  X  10  “ 17  T102  exp(203.5/RT) 
cm3/(molecule-s) 
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COMMENTS 


Jurgen  Troe,  University  of  Gottingen,  Germany.  I  would 
like  to  remind  you  that  the  reaction  C2H5  +  02,  which 
you  describe,  corresponds  completely  to  the  reaction  CO 
+  OH,  that  is,  biexponential  kinetics,  pressure-dependent, 
pseudo-second-order  rate  coefficients,  and  low  and  high 
pressure  limiting-ranges  are  observed.  In  my  analysis  of  the 
CO  +  OH  reaction  [1],  I  came  to  the  conclusion  that  a 
two-dimensional  treatment  with  separate  energy  (E)  and 
angular  momentum  (J)  dependencies  of  the  rate  constants 
is  required,  in  particular  if  a  loose  entrance/rigid  exit  com¬ 
plex-forming  birnolecular  reaction  is  considered.  Can  your 
code  handle  this?  Another  question:  what  energy  transfer 
parameters  did  you  use? 
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Author’s  Reply.  The  treatment  you  use  in  your  paper  is 
equivalent  to  what  we  call  pVT-J  theory  (microcanonical 
variational  theory  with  J  conservation),  with  a  pseudo¬ 
strong-collider  approximation  for  stabilization.  We  cer¬ 
tainly  do  have  such  a  capability.  We  conclude  that  J  con¬ 
servation  is  not  a  very  important  factor  in  the  present 
reaction  (see  Fig.  5  of  the  paper).  Our  full  one-dimensional 
master  equation  treatment  goes  beyond  such  a  simple 
treatment,  at  least  as  long  as  J  conservation  is  not  a  major 
concern. 

For  lack  of  anything  better  to  use  at  this  point,  we  as¬ 
sumed  a  single-exponential  down  model  for  the  energy 


transfer  function  with  /il\\  =  200  cm  .  We  arrived  at  this 
number  by  testing  against  some  thermally  activated  disso¬ 
ciation/recombination  reactions  with  helium  as  the  bath 
gas.  It  worked  well  in  the  present  calculations,  so  we  did 
not  alter  it. 


Horst  Hippier,  University  of  Karlsruhe,  Germany.  Col- 
lisional  energy  transfer  controls  the  transition  between 
pressure-independent  low-  and  high-pressure  limiting  rate 
constants.  How  has  this  been  treated,  and  how  sensitive  is 
the  product  distribution  on  the  collisional  energy  transfer 
parameters  in  this  range?  Secondly,  at  high  temperatures 
the  kinetic  system  converts  into  a  reaction  out  of  a  pre¬ 
equilibrium.  Therefore,  the  rate  constant  for  H02  forma¬ 
tion  is  not  only  controlled  by  the  properties  of  TS2  but  also 
by  the  equilibrium  constant  for  C2H5  +  02  <-»  C2H502. 
Thus,  also  thermochemistry  of  C2H502  equilibrium  con¬ 
trols  the  reaction  rate. 

Author’s  Reply.  Regarding  the  first  part  of  your  question, 
see  the  response  to  the  question  by  Jurgen  Troe.  Secondly, 
I  do  not  know  what  you  mean  by  “preequilibrium”.  As  al¬ 
luded  to  in  the  paper,  the  rate  coefficient  in  the  high-tem- 
perature  regime  does  not  depend  on  the  thermochemistry 
of  C2H502.  However,  the  temperature  at  which  the  high- 
temperature  regime  is  reached  does.  Specifically,  the  tem¬ 
perature  at  which  the  transition  regime  occurs  depends 
significantly  on  the  C2H502  bond  energy. 


